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A  two-dimensional  two-phase  mass  transport  model  has  been  developed  to  predict  methanol  and  water 
crossover  in  a  semi-passive  direct  methanol  fuel  cell  with  an  air-breathing  cathode.  The  mass  transport 
in  the  catalyst  layer  and  the  discontinuity  in  liquid  saturation  at  the  interface  between  the  diffusion  layer 
and  catalyst  layer  are  particularly  considered.  The  modeling  results  agree  well  with  the  experimental  data 
of  a  home-assembled  cell.  Further  studies  on  the  typical  two-phase  flow  and  mass  transport  distributions 
including  species,  pressure  and  liquid  saturation  in  the  membrane  electrode  assembly  are  investigated. 
Finally,  the  methanol  crossover  flux,  the  net  water  transport  coefficient,  the  water  crossover  flux,  and  the 
total  water  flux  at  the  cathode  as  well  as  their  contributors  are  predicted  with  the  present  model.  The 
numerical  results  indicate  that  diffusion  predominates  the  methanol  crossover  at  low  current  densities, 
while  electro-osmosis  is  the  dominator  at  high  current  densities.  The  total  water  flux  at  the  cathode  is 
originated  primarily  from  the  water  generated  by  the  oxidation  reaction  of  the  permeated  methanol  at 
low  current  densities,  while  the  water  crossover  flux  is  the  main  source  of  the  total  water  flux  at  high 
current  densities. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cells  (DMFCs)  have  been  considered  as  one 
of  the  most  promising  power  sources  for  portable  devices  due  to 
its  high  energy  density,  simple  system,  fast  refueling,  low  operat¬ 
ing  temperature,  and  system  simplicity.  Fuel  cells  with  neither  fuel 
pumps  nor  oxidant  suppliers  are  regarded  as  ‘passive’,  while  those 
with  only  one  reactant  by  forced  convection  as  ‘semi-passive’.  In 
particular,  air-breathing  DMFCs,  which  are  exposed  to  atmosphere 
at  cathode  side  and  draw  air  for  its  operation  entirely  through  natu¬ 
ral  convection,  have  attracted  much  attention  since  simpler  system 
eliminates  the  parasitic  power  losses  caused  by  auxiliary  devices 
[1-8].  However,  the  wide  application  of  air-breathing  DMFCs  is 
still  hindered  by  many  practical  problems  in  materials,  fabrica¬ 
tion,  and  operation.  Methanol  crossover  and  water  crossover  are 
considered  as  two  of  the  major  problems  in  the  DMFCs  opera¬ 
tion.  Methanol  crossover  will  increase  the  cathode  polarization 
and  lower  both  cell  performance  and  energy  conversion  efficiency. 
Water  management  will  cause  flooding  at  the  cathode  due  to  dif¬ 
ficulties  in  removing  excessive  water  and  further  degrade  the  cell 
performance.  The  understanding  on  the  mechanisms  and  contribu¬ 
tors  to  the  methanol  crossover  and  the  water  crossover  will  promote 
the  application  of  the  air-breathing  DMFCs. 
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Recently,  many  experimental  investigations  have  been  carried 
out  for  passive  and  air-breathing  DMFCs  on  aspects  including 
optimizing  operation  conditions  [2-8],  new  fuel  delivery  systems 
[9-11],  new  structures  of  membrane  electrode  assembly  (MEA) 
[12,13],  effect  of  methanol  crossover  [14],  and  water  manage¬ 
ment  [15,16].  However,  there  are  still  many  difficulties  to  perform 
detailed  in  situ  measurement  in  an  operating  cell.  The  numeri¬ 
cal  simulation  provides  an  available  tool  to  deeply  understand  the 
transport  interaction  and  electrochemical  phenomena  inside  the 
air-breathing  DMFC. 

Numerous  mathematical  models  have  been  reported  for  air- 
breathing  PEMFCs  [17-27]  and  DMFCs  [18-34],  among  which,  most 
considered  only  single  phase  model  [17-30],  while  a  few  adopted 
two-phase  and  multi-phase  approaches  [31-34].  Some  papers 
modeled  heat  transport  in  the  fuel  cells  [18,21-27,29,31,33,34], 
though  other  papers  neglected  the  thermal  effect.  Most  of  the 
models  studied  the  vertical  fuel  cells  except  Xiao  and  Faghri  [34] 
in  which  the  cell  was  orientated  horizontally.  In  addition,  some 
of  the  previous  studies  treated  the  catalyst  layer  as  an  interface 
[21,23,27-29,31],  others  considered  the  mass  transport  in  the  cat¬ 
alyst  layer.  In  these  models,  the  liquid  saturation  was  simulated  as 
continuous  at  the  interface  between  the  diffusion  layer  and  catalyst 
layer.  Actually,  a  discontinuity  in  liquid  saturation  exists  at  the  inter¬ 
face  between  two  kinds  of  porous  media  due  to  continuous  capillary 
pressure  and  different  physical  properties  according  to  the  classi¬ 
cal  multi-phase  flow  theory  in  porous  media.  Recently,  only  Xu  et 
al.  [35]  and  Pasaogullari  and  Wang  [36]  accounted  for  the  sudden 
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Nomenclature 

Av  specific  area  (m2  m-3 ) 

C  molar  concentration  (mol  m-3 ) 

C0  distribution  parameter 

D  diffusivity  (m2  s-1 ) 

/  friction  factor 

F  Faraday  constant  (C  mol-1 ) 

g  gravitational  acceleration  (ms-2) 

ha  height  of  anode  flow  channel  (m) 

hm  mass  transfer  coefficient  (ms-1) 

/  local  current  density  (A  m-2 ) 

JCell  cell  current  density  (A  m-2 ) 

IP  parasitic  current  density  (A  m-2 ) 

ja  anode  current  density  (A  m-3 ) 

jc  cathode  current  density  (A  m-3 ) 

jo  exchange  current  density  (A  m-2 ) 

kc  condensation  rate  (s-1 ) 

/<e  evaporation  rate  (atm-1  s-1 ) 

kr  relative  permeability 

I<  permeability  of  porous  media  (m2 ) 

l  length  of  flow  channel  (m) 

m  mass  flux  (kg  m-2  s-1) 

m  source  term  in  mass  conservation  equation 

(kgm-3  s-1) 

M  molecular  weight  (kg  mol-1 ) 

nd  electro-osmotic  drag  coefficient 

N  mol  flux  (mol  m-2  s-1),  or  Nation 

p  pressure  (Pa) 

pc  capillary  pressure  (Pa) 

pg  gas  phase  pressure  (Pa) 

Pi  liquid  phase  pressure  (Pa) 

R  gas  constant  (J  mol-1 1<-1 ) 

RH  relative  humidity 

R  source  term  in  species  conservation  equation 

(molm-3  s-1) 

^contact  ohmic  contact  resistance  ( £2  m2 ) 

Rw  interfacial  transfer  rate  of  water  between  liquid  and 
vapor  (molm-3  s-1) 
s  liquid  saturation 

Sc  Schmidt  number 

Sh  Sherwood  number 

T  temperature  (I<) 

u  superficial  velocity  vector  (ms-1) 

u  phase  velocity  (ms-1) 

Ugj  drift  velocity  (ms-1) 

Veen  cell  voltage  (V) 

Vo  thermodynamic  equilibrium  potential  (V) 

wa  width  of  anode  flow  channel  (m) 

x  coordinate  (mm),  or  mole  fraction  in  liquid  solution 

y  coordinate  (mm),  or  mole  fraction  in  gas  mixture 

Greek  letters 

a  void  fraction  in  flow  channel 

aa  transfer  coefficient  of  anode 

ac  transfer  coefficient  of  cathode 

aw  net  water  transport  coefficient 

y  reaction  order 

8  thickness  of  porous  layer  (m) 

s  porosity  of  porous  layer 

p  overpotential  (V) 

0C  contact  angle  (°) 

k  ionic  conductivity  of  membrane  (£2-1  m-1 ) 

fi  viscosity  (kg  m-1  s-1 ) 


p  density  (kgm-3) 

a  surface  tension  (Nm-1) 

Superscripts 

acl  anode  catalyst  layer 

adl  anode  diffusion  layer 

afc  anode  flow  channel 

amb  ambient 

ccl  cathode  catalyst  layer 

cdl  cathode  diffusion  layer 

eff  effective  value 

mem  membrane 

ref  reference  value 

sat  saturated  value 

—  average  value 

Subscripts 
a  anode 

c  cathode 

Cell  cell 

g  gas  phase 

1  liquid  phase 

MeOH  methanol 

02  oxygen 

v  vapor 

w  liquid  water 


change  of  liquid  saturation  in  their  models  on  the  water  transport 
in  DMFCs  and  PEMFCs,  respectively. 

In  the  present  study,  a  two-dimensional  two-phase  mass  trans¬ 
port  model  is  developed  to  predict  methanol  and  water  crossover 
in  a  semi-passive  direct  methanol  fuel  cell  with  an  air-breathing 
cathode.  In  addition  to  the  mass  transport  in  the  catalyst  layer, 
the  discontinuity  in  liquid  saturation  at  the  interface  between  the 
diffusion  layer  and  catalyst  layer  is  also  taken  into  account.  The 
drift-flux  model  is  adopted  to  describe  the  two-phase  flow  in  the 
anode  channels.  The  mass  transport  in  the  anode  porous  region  is 
modeled  by  the  classical  multi-phase  flow  in  porous  media  theory, 
in  which  various  phases  are  regarded  as  distinct  fluids  with  indi¬ 
vidual  thermodynamic  and  transport  properties  and  with  different 
flow  velocities.  Flowever,  the  mass  transport  in  the  cathode  is  for¬ 
mulated  based  on  the  unsaturated  flow  theory  (UFT),  in  which  the 
gas  pressure  is  assumed  to  be  constant  across  the  entire  porous 
region.  The  predicted  results  are  validated  with  experimental  data 
by  a  home-assembled  cell.  Typical  two-phase  flow  and  mass  trans¬ 
port  distributions  at  the  anode  and  the  cathode,  methanol  crossover 
flux,  net  water  transport  coefficient,  water  crossover  flux,  as  well 
as  the  total  water  flux  at  the  cathode  and  their  contributors  are 
investigated  with  this  model. 

2.  Numerical  model 

A  vertical  air-breathing  DMFC  simulated  in  the  present  study  is 
shown  in  Fig.  1  and  methanol  aqueous  solution  is  fed  by  a  pump 
in  the  direction  as  marked.  The  computational  domain  consists  of 
all  six  parts  of  the  air-breathing  DMFC,  including  anode  flow  chan¬ 
nel  (AFC),  anode  diffusion  layer  (ADL),  anode  catalyst  layer  (ACL), 
membrane  (MEM),  cathode  catalyst  layer  (CCL)  and  cathode  dif¬ 
fusion  layer  (CDL).  At  the  anode,  methanol  aqueous  solution  is  fed 
into  the  AFC  and  then  transfers  through  the  ADL  into  the  ACL,  where 
one  part  of  methanol  is  oxidized  to  generate  C02  gas,  while  another 
part  permeates  through  the  MEM  into  the  cathode.  The  produced 
C02  gas  goes  through  the  ADL,  emerges  into  the  AFC  in  bubbles  and 
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Fig.  1.  Schematic  of  a  vertical  air-breathing  DMFC. 


then  removes  with  flowing  methanol  solution  in  the  AFC.  At  the 
cathode,  the  oxygen  from  the  ambient  air  diffuses  through  the  CDL 
into  the  CCL,  where  it  reduces  with  the  electrons  from  the  circuit 
and  the  protons  from  the  anode  to  form  water.  The  produced  water 
is  expelled  through  the  CDL  by  the  capillary  pressure  along  with 
evaporation. 

Several  coupled  complex  physicochemical  processes  occur 
simultaneously  in  such  an  operation  air-breathing  DMFC,  which 
makes  it  quite  difficult  to  simulate  comprehensively.  For  the  numer¬ 
ical  realization  of  the  present  model,  the  following  assumptions, 
which  are  reasonable  and  acceptable  to  some  extent  based  on  the 
practical  operation  DMFC,  are  used: 

(1 )  The  cell  operates  under  steady-state  and  isothermal  conditions. 

(2)  Thermal  physical  properties  of  fluids  are  considered  as  constant. 

(3)  The  fluid  flow  in  the  AFC  is  one-dimensional,  i.e.,  the  flow  along 
x-direction  is  ignored. 

(4)  Considering  the  thicknesses  of  ADL,  ACL,  MEM,  CCL  and  CDL  are 
much  smaller  than  their  sizes  along  y-direction,  the  flow  and 
mass  transport  along  y-direction  are  neglected. 

(5)  The  ADL,  ACL,  MEM,  CCL  and  CDL  are  isotropic  porous  media. 

(6)  At  the  anode,  the  gas  fluid  is  considered  only  from  the  produced 
CO2  gas,  i.e.,  the  evaporation  of  methanol  and  water  is  ignored. 

(7)  The  methanol  permeated  from  the  anode  is  oxidized  completely 
at  the  interface  between  the  MEM  and  CCL. 

(8)  The  gas  pressure  at  the  cathode  is  a  constant  since  the  cathode 
is  exposed  to  the  ambient. 

(9)  At  the  cathode,  the  gases  are  transferred  from  the  ambient  to 
the  CDL  by  natural  convection. 

2.1.  Anode  flow  channel 

Visualization  experiments  for  DMFCs  [37,38]  showed  that  the 
two-phase  flow  pattern  in  AFC  varied  from  homogeneous  to  slug 
flow  depending  on  the  cell  operating  conditions.  The  amount  of  gas 
C02  in  the  AFC  increased  with  increasing  current  density,  and  the 
gas  bubbles  or  slugs  were  removed  periodically  by  the  solution  in 
the  channel.  Thus  it  can  be  seen  that  a  complicated  two-phase  flow 
happens  in  the  AFC  of  a  practical  operating  DMFC.  In  the  present 
work,  a  drift-flux  model,  which  is  firstly  used  for  a  DMFC  by  Wang 


and  Wang  [39],  is  adopted  to  describe  gas-liquid  two-phase  flow  in 
the  anode  channel.  Details  are  presented  below. 

The  continuity  equations  for  two  phases  in  the  AFC  can  be  given 
by: 

mfc 

V[aui(1  -  a)]  =  — (liquid  phase),  (1) 

mafc 

V(pgUgO!)= — f-  (gas  phase),  (2) 

where  the  subscripts  T  and  ‘g’  refer  to  the  liquid  phase  and  the 
gas  phase,  respectively,  u  represents  the  average  velocity  across  the 
flow  channel,  a  is  the  cross-sectional  mean  gas  void  fraction,  mafc 
denotes  the  mass  flux  at  the  interface  between  the  AFC  and  ADL, 
and  ha  is  the  height  of  the  AFC  cross-section. 

The  momentum  equation  for  the  vertical  anode  channel  can  be 
written  as: 

f 

VtouiiiKl  -  a)]  +  V(pgUgUga)  =  -Vp  -  2 pm^-umum  -  pmg,  (3) 

where  p  is  the  liquid  pressure  in  the  AFC,  pm  the  two-phase  volume- 
averaged  density  pm  =  p\{  1  -a)  +  pga,fm  the  friction  factor,  Dh  the 
hydraulic  diameter  of  the  AFC,  and  um  represents  the  average  veloc¬ 
ity  of  the  two-phase  mixture  and  is  given  as: 


um  =  aug  +  (  1  -a)u\.  (4) 

According  to  the  drift-flux  model,  the  relationship  between  the 
velocities  of  liquid  and  gas  phases  can  be  described  by  the  following 
equation: 


ug  =  Co[aug  +  (1  -  a)U\]  +  ugj,  (5) 

where  Co  represents  the  distribution  parameter  and  is  the  drift 
velocity.  The  distribution  parameter  and  the  drift  velocity  for  slug 
flow  through  a  rectangular  channel  are  given  by: 


C0  =  1 .35  —  0.35,  /  —  , 
Pi 


«a.(°.23  +  °.,3  £)/ 


(A  ~  Pg)gWa 

Pi 


(6) 

(7) 


where  wa  is  the  width  of  the  AFC. 

The  species  balance  equation  for  methanol  in  the  AFC  can  be 
written  as: 

Nafc 

V[(l-a)u,CMeoH]  = - (8) 


where  CMeoH  is  the  cross-sectional  average  methanol  concentration 
and  N^oh  represents  the  mole  flux  of  methanol  at  the  interface 
between  the  AFC  and  ADL. 


2.2.  Anode  porous  region 

The  anode  porous  region  consists  of  the  ADL  and  ACL.  The 
ADL  is  composed  of  inactive  carbon  and  no  electrochemical  reac¬ 
tion  occurs  within  this  region.  The  ACL  provides  catalytic  sites  for 
methanol  oxidation.  The  transport  behavior  in  the  anode  porous 
region,  where  the  transport  of  methanol  solution  is  driven  by  con¬ 
vection  and  diffusion,  is  modeled  based  on  the  classical  two-phase 
flow  in  porous  media.  The  mass  conservation  equations  for  the  two 
phases  in  the  anode  porous  region  are  given  by: 

V(auD  =  mu  (liquid  phase),  (9) 

V(pgUg)  =  rhga  (gas  phase),  (10) 

where  u  is  the  superficial  velocity  vector  in  the  porous  media  and 
rh  represents  the  mass  source  (see  Table  1),  subscript  ‘a’  refers  to 
the  anode. 


Table  1 

Correlations  used  in  this  model. 
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Parameters 


Expressions 


Relative  permeability 

Effective  diffusion  coefficients  of  species 


General  generation  rate  of  mass  in  liquid  phase 


General  generation  rate  of  mass  gas  phase 

Mole  generation  rate  of  species 

The  saturation  pressure  of  water  vapor 

Interfacial  transfer  rate  of  water  between  liquid  and  vapor 


kr\  =  s3  liquid 
fcrg  =  (l-s)3  gas 

Dfff  =  D,e1-5s1-5  i:  02,  vapor  CDL/CCL 

^MeOH.lS1’5^1’5 

Deff  - 

MeOH 


S  +  e  N 


mu  = 


mlc  = 


™g.a  =  ■ 


ADL 

ACL 


[^/(DMeOH.l^1-5^1-5)  +  &N  /  (f^MeOH,N  1  5  )  ] 

f^MeOH.N^1’5  MEM 

'  0  ADL 

-MwjL -MMe0Hj|=  ACL 


-MwRw 


CDL 


,Mw(^-6^)  -Mw*w  CCL 

r  0  ADL 

ACL 

r  0  ADL 


^MeOH  =  S  /’ 

\~W  ACL 


Rq2  = 


~  4  F 


CDL 

CCL 


logioPvat  =  -2.1794  +  0.02953 (T  -  273)  -  9.1837  x  10"5(T  -  273  f  +  1.4454  x  1CT7(T  -  273 )3  atm 
yvpg  <  Pvat 


Rw  — 


k^y,P%-P?) 

_  kc  +L_£+(yvpg  _ psa«)  yvPg  >  psa, 


The  momentum  equations  in  the  anode  porous  region  are 
described  by  Darcy’s  law  as: 

m  =  -/C  — Vpla  (liquid  phase),  (11) 

Mi 

k 

ug  =  -/C^Vpg,a  (gas  phase),  (12) 

Mg 

where  K,  and  krg  denote  the  absolute  permeability  of  porous 
media,  the  relative  permeability  of  liquid  and  gas  (see  Table  1), 
respectively,  and  /xg  are  the  viscosity  of  liquid  and  gas.  The 
gas  pressure,  pg,  and  the  liquid  pressure,  relate  to  the  capillary 
pressure,  pc,  which  is  defined  by: 

/£n!/2 

Pc  =  Pg  -  Pi  =  &  cos  0cl- )  ■/(+  (13) 

where  a  is  the  surface  tension  of  methanol  solution,  0C  the  contact 
angle,  e  the  porosity  of  porous  media,  and  J(s)  stands  for  the  Levertte 
function  and  can  be  given  by  [40]: 


2.3.  Membrane 

Only  mass  transport  of  the  liquid  methanol  and  water  across  the 
MEM  is  taken  into  account  due  to  the  impermeability  of  the  proton 
exchange  membrane  for  gases.  It  is  noted  that  the  methanol  trans¬ 
port  through  the  membrane  is  driven  by  diffusion,  electro-osmosis, 
and  hydraulic  pressure  gradient  across  the  membrane.  Thus,  the 
flux  of  methanol  crossover  can  be  given  by: 

\jmem  _  nmem-effv7/^mem  ,  „  _  ^  /-mem  ^XPl.c  —  Pi, a) 

iVMeOH“  ^MeOH  v L'MeOH‘1"“d,MeOH  c  LMeOH  ..  s  ’  ^  1DJ 

r  /xiOmem 

where  D^Q^eff  is  the  effective  diffusivity  of  methanol  in  the  mem¬ 
brane,  (Wm  the  thickness  of  the  membrane,  ndiMe0H  denotes  the 
electro-osmotic  drag  coefficient  of  methanol,  nd)Me0H  =  nd,w*MeOH» 
nd  w  the  electro-osmotic  drag  coefficient  of  water,  and  xMeoH  the 
mole  fraction  of  methanol  solution. 

The  methanol  flux  is  uniform  through  the  membrane  due  to  no 
electrochemical  reaction  in  the  MEM,  i.e., 


,  f  1.417(1  -s)-  2.120(1  -s)2  +  1.263(1  -s)3  0  <  0C  <  90°  _ 

J(s)  =  s  (14) 

1  1.417s  -  2.120s2  +  1.263s3  90°  <  0C  <  180° 

The  wettability  of  the  ADL  and  ACL  is  set  to  be  hydrophobic  in 
the  present  study,  hence  the  contact  angles  larger  than  90°  (see 
Table  4). 

The  species  conservation  of  methanol  in  the  anode  porous  region 
can  be  expressed  as: 

VtUjCjvieOH)  +  V(-D^OHVCMeOH)  =  ^MeOH  ?  (15) 


VN{JeeS1H  =  0.  (17) 

Furthermore,  the  transport  mechanism  of  liquid  water  through 
the  membrane  is  similar  to  that  of  methanol.  However,  the  water 
concentration  difference  across  the  membrane  is  considerably 
small  leading  to  the  water  diffusion  being  neglected.  Therefore,  the 
water  flux  is  expressed  as  follows: 


jumem 

iVw 


A  K(Pl,c-Pl,a) 
Mw  Ml^mem 


(18) 


where  Cmcoh  and  D^0H  represent  the  methanol  concentration 
and  the  effective  diffusion  coefficient  of  methanol  (see  Table  1), 
respectively.  The  first  term  on  the  left-hand  side  of  Eq.  (15)  denotes 
the  convection  term  of  methanol  transport,  while  the  second  term 
describes  the  diffusive  transport,  and  the  source  term  ^MeOH  stands 
for  the  mole  generation  rate  of  methanol  due  to  the  electrochemical 
reaction  (see  Table  1 ). 


2.4.  Cathode  porous  region 

The  cathode  porous  region  consists  of  the  CDL  and  CCL.  Simi¬ 
larly,  the  CDL  is  composed  of  inactive  carbon  and  no  electrochemical 
reaction  occurs  within  this  region,  while  the  CCL  provides  catalytic 
sites  for  oxygen  reduction.  The  wettability  of  the  CDL  and  CCL  is  also 
set  to  be  hydrophobic.  The  two-phase  flow  in  the  cathode  porous 
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region  is  modeled  with  the  unsaturated  flow  theory.  Since  gas  pres¬ 
sure  at  the  cathode  is  assumed  to  be  a  constant,  the  transport  of 
gases  in  the  cathode  porous  region  is  dominated  by  diffusion.  The 
liquid  water  generated  in  the  CCL  removes  through  the  CDL  by  the 
capillary  pressure  coupled  with  evaporation. 

The  mass  conservation  of  the  liquid  in  the  cathode  porous  region 
can  be  modeled  by: 

V(pwu,,c)  =  m,.c;  (19) 

where  subscript  ‘w’  refers  to  liquid  water  and  ‘c’  the  cathode. 

The  momentum  equation  of  the  liquid  can  be  obtained  from 
Darcy’s  law: 

«l,c  =  -^^£Vp1,c>  (20) 

where  /xw  is  the  viscosity  of  the  liquid  water.  With  the  aid  of 
the  definition  of  the  capillary  pressure  shown  in  Eq.  (13)  and  the 
assumption  of  the  constant  gas  pressure  across  the  cathode  porous 
region,  Eq.  (20)  can  be  rewritten  as: 

uliC  =  K^awcos0c(|)1/2VJ(s).  (21) 

The  mass  transport  of  the  gases  in  the  cathode  porous  region 
can  be  modeled  by: 

V(-Df  VCv)  =  Rw,  (22) 

V(-DefVCo2)  =  Ro2,  (23) 

where  DJff  and  represent  the  effective  diffusivity  of  vapor  and 

oxygen,  respectively,  Rw  is  the  interfacial  transfer  rate  of  water 
between  liquid  and  vapor,  and  Rq2  the  mole  generation  rate  of 
oxygen  (see  Table  1). 

2.5.  Electrochemical  kinetics 


At  the  anode,  a  simplified  Tafel  equation  is  employed  to  model 
the  kinetics  of  methanol  oxidation  reaction  (MOR): 


where  Av,a  represents  the  specific  area  of  the  anode,  j’oeMe0H  *s  the 
reference  exchange  current  density  of  the  anode,  CJ^0H  the  refer¬ 
ence  methanol  concentration,  and  y  the  reaction  order  and  given 
by: 


(25) 


which  means  the  methanol  oxidation  is  a  zero-order  reaction  when 
methanol  concentration  is  higher  than  a  reference  value,  otherwise, 
a  first-order  reaction  is  considered  in  this  model. 

At  the  cathode,  the  first-order  Tafel  kinetics  is  used  to  describe 
the  oxygen  reduction  reaction  (ORR): 


jc  —  (1  —  s)A/,doeo2 


(26) 


where  AV)C  represents  the  specific  area  of  the  cathode,  the  ref¬ 
erence  exchange  current  density  of  oxygen,  and  C^f  the  reference 
oxygen  concentration. 


2.6.  Boundary  conditions 


At  the  inlet  of  the  AFC,  all  variables  can  be  given  according  to  the 
fuel  supply  conditions: 


r,afc 


nflaWa  ’ 


^afc  _  f'in 
LMeOH  “  LMeOH’ 


Pi=PiS 


(1  -  a)  =  s  =  sjna. 

(27) 


The  continuity  and  mass/species  flux  balance  conditions  are 
employed  at  each  interface.  Detailed  description  of  the  boundary 
conditions  for  each  interface  is  given  in  Table  2.  It  should  be  noted 
that  the  capillary  pressure  is  continuous  at  x  =  xadi  and  x  =  xcc\ ,  while 
different  properties  of  the  porous  media  result  in  a  discontinuity  of 
the  liquid  saturation  at  these  interfaces. 

It  is  difficult  to  determine  the  liquid  saturation  at  the  outer  sur¬ 
face  of  the  CDL  in  the  air-breathing  DMFC.  Actually,  it  changes  with 
the  operation  condition  and  the  current  density.  For  simplification, 
a  constant  value  is  given  to  the  liquid  saturation  at  the  CDL  surface 
in  present  model,  i.e.,  s°c  =  0.05  [31  ]. 


Table  2 

Boundary  conditions. 


Interface 

Species 

Flux 

Pressure 

Liquid  saturation 

x  =  0 

f'Ofc  r'adl 

MeOH  ~  ^ MeOH 

juafc  _  |Vadl 

MeOH  —  MeOH 

rnafc  —  madl 
ml  "‘MeOH 

mafc  —  madl 

"g  "‘MeOH 

pfc=pf' 

(1  -  a)  =  sfc  =  sfdl 

*  =  *adl 

/-adl  /-acl 

'“MeOH  —  '“MeOH 

Madl  _  ivracl 

‘*MeOH  —  ‘*MeOH 

o. 

II 

J3 

cos  0adl(jg)V2/(sfdl: 

l  =  C°S  <?ad(g 

\  1/2 

)  J(sf) 

madl  =  macl 

II 

madl  =  macl 

*  =  *acl 

/-acl  _  /-mem 

MeOH  '“MeOH 

/yacl  _  lumem 

MeOH  MeOH 

vp?ci  =  -^«) 

m1ad  =  m[nem,  macl  =  0 

o 

II 

a? 

> 

^  =  ^mem 

/-mem  /-ccl  n 

MeOH  —  '“MeOH  —  u 

ivrmem  _  iwccl 

MeOH  —  MeOH 

Vpfcl|-w(mfcl) 

Ng  =  0,  =  0,  m™em  =  mfcl,  rnf  =  0 

o 

II 

> 

X  =  Xcd 

CP 

II 

CP 

mccI  —  Mcdl 

JV02  “  02 

s. 

II 

^3 

cos  Pccl(|i)1/2j(sfd) 

=  COSPcdl(^) 

i’V) 

q 

ii 

Ncd  =  Ncdl 

mcd  =  mf 
mf  =  mf 

2: 

II 

^3 

X  =  Xcdl 

Ng  =  ho2(C%*-Cg) 

N'f  =  hv(Qmb  -  Cvcdl) 

Pgdl=Pe,c 

<~cdl  <-0 

_  \c 
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The  gases  transfer  between  the  CDL  and  the  ambient  air  by 
natural  convection.  Therefore  the  mass  transfer  coefficient  at  the 
interface  of  x  =  xcdI  can  be  obtained  from  the  natural  convection 
correlation  on  a  vertical  surface  [17]: 


Sh=^=  0.68  + 


0.67  Ra, 


1/4 

M 


[1  +  (0.492/Sc )' 


,9/16l4/9’ 


(28) 


where  Sh  is  the  Sherwood  number,  Sc  the  Schmidt  number  (Sc  =  p/D), 
and  Rqm  the  Rayleigh  number  (Rcim  =  (yg(ww  -  wOQ)L3)lvD). 


2.7.  Cell  performance 


The  electrochemical  reaction  rate  varies  with  the  reactant  con¬ 
centration  along  both  x-direction  and  y-direction.  As  a  result,  the 
local  current  density  and  the  average  cell  current  density  can  be 
given  by: 


=  /  Jadx, 

J  acl 

J'/ac^dxdy 


1  = 
and 


/ceil  = 


/ 


(29) 


(30) 


Furthermore,  the  parasitic  current  density  caused  by  the  methanol 
crossover  is  calculated  by: 


Ip  =  6FNJJJ|qH. 

The  current  balance  in  the  cell  is  described  as: 


I  +  Ip  =  /  jcdx. 

J  ccl 

Thus,  the  cell  voltage  is  expressed  as: 


''Cell 


—  Vq  —  Pa  —  Pc  —  /cell  (  /^contact  H - — — 


(31) 

(32) 

(33) 


where  Vo  is  the  thermodynamic  equilibrium  potential  of  a  DMFC, 
^contact  the  contact  resistance  of  the  fuel  cell,  and  /c  denotes  the 
proton  conductivity  of  the  membrane. 


2.8.  Numerical  method 


Table  3 

Operating  parameters. 


Parameter 

Symbol 

Value 

Unit 

Anode  flow  channel  length 

l 

2.1  x  10“2 

m 

Anode  flow  channel  height 

K 

2.0  x  10-3 

m 

Anode  flow  channel  width 

wa 

1.8  x  10“3 

m 

Anode  flow  channel  number 

n 

7 

- 

Thickness  of  ADL  and  CDL 

Jdl 

2.6  x  10-4 

m 

Thickness  of  ACL  and  CCL 

Cl 

2.0  x  10-5 

m 

Thickness  of  MEM 

lm 

1.8  x  10“4 

m 

Operating  temperature 

T 

298.0 

I< 

Anode  inlet  pressure 

Pt 

1.013  x  105 

Pa 

Cathode  gas  pressure 

Pg,C 

1.013  x  105 

Pa 

Total  anode  flow  rate 

0 

1.0 

ml  min-1 

Inlet  methanol  concentration  at  anode 

/-in 

'“'MeOH 

4000 

mol  m3 

Relative  humidity  of  air 

RHair 

60% 

- 

Contact  resistance 

^contact 

o 

X 

o 

00 

Q  m2 

anode  and  the  cathode.  The  anode  catalyst  loading  was  3.0  mg  cm-2 
with  Pt-Ru  black  (1:1,  a/o),  while  the  cathode  catalyst  loading  was 
5.0  mg  cm-2  using  Pt  black.  The  membrane  electrode  assembly  was 
formed  by  hot  pressing  a  pretreated  Nation  117  membrane  between 
the  two  electrodes  at  135  °C  and  10  MPa  for  3  min.  The  prepared 
MEA  was  sandwiched  between  two  stainless  steel  plates,  in  which 
the  straight  flow  channels  were  fabricated.  The  methanol  solution 
was  fed  into  the  anode  by  a  peristaltic  pump.  The  cathode  electrode 
was  exposed  to  atmosphere  and  oxygen  was  drawn  directly  from 
the  ambient  air  by  natural  convection.  The  voltage-current  curves 
were  recorded  by  an  Arbin  Fuel  Cell  Test  system  (FCTs-1000). 

The  modeling  results  of  the  cell  polarization  curve  are  compared 
with  the  experimental  data  in  Fig.  2.  It  can  be  seen  that  the  numer¬ 
ical  results  agree  well  with  the  experimental  data  within  a  wide 
current  density  range.  The  slight  disagreement  at  low  current  den¬ 
sities  might  be  resulted  from  the  Tafel  kinetics  employed  in  this 
model.  Based  on  the  good  agreement  with  the  data,  the  further 
studies  on  species  distributions,  pressure,  liquid  saturation,  as  well 
as  methanol  and  water  crossover  are  performed  as  follows  with  the 
help  of  the  above  model.  Unless  specified  otherwise,  the  modeling 
results  on  the  distributions  are  obtained  with  4  M  methanol  con¬ 
centration  and  at  the  current  density  of  720.78  A  rrr2  ( rja  =  0.385  V). 


The  above  governing  equations  are  iteratively  solved  by  the  finite 
volume  method  using  a  self-written  code  in  FORTRAN  language 
based  on  the  tridiagonal  matrix  algorithm  (TDMA)  until  the  con¬ 
vergence  criterions  are  satisfied.  The  operating  parameters  and  the 
physicochemical  properties  used  in  the  calculations  are  listed  in 
Tables  3  and  4,  respectively.  For  a  given  anode  overpotential 
the  methanol  concentration  distribution  can  be  obtained  using  Eqs. 
(1 )— (18).  Then  the  average  cell  current  density,  /Cen,  and  the  para¬ 
sitic  current  density,  /p,  can  be  determined  from  Eqs.  (30)  and  (31). 
The  cathode  overpotential,  pc,  can  be  calculated  from  Eq.  (32)  based 
on  the  oxygen  concentration  distribution  and  the  cathode  current 
density,  jc,  obtained  from  Eqs.  (19)-(23)  and  (26).  Finally,  the  anode 
and  the  cathode  overpotentials  are  substituted  into  Eq.  (33)  to  cal¬ 
culate  the  cell  voltage  VCell.  Such  a  computation  process  is  repeated 
until  the  polarization  curve  is  obtained. 

3.  Results  and  discussion 

3.1.  Model  validation 

An  air-breathing  DMFC  with  an  active  area  of  2.1  x  2.1  cm2  was 
fabricated  and  tested  to  validate  the  above  model.  A  carbon  paper 
(Toray,  TGP-H-090)  with  15%  polytetrafluoroethylene  (PTFE)  wet¬ 
proofing  treatment  was  used  as  the  backing  support  layer  in  the 


3.2.  Reactants  concentration  distributions 

It  is  known  that  methanol  transfers  from  the  AFC  through  the 
ADL  to  the  ACL  at  the  anode  and  permeates  through  the  MEM  to 
the  cathode.  Fig.  3  shows  the  methanol  concentration  distribution 
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Table  4 

Physicochemical  properties. 


in  the  anode  porous  region  and  MEM.  Apparently,  the  methanol 
concentration  decreases  linearly  across  the  ADL,  while  a  non-linear 
drop  is  observed  across  the  ACL  due  to  the  methanol  oxidation 
reaction  there.  Finally,  it  goes  down  to  zero  at  x  =  xmem  since  the  per¬ 
meated  methanol  through  the  MEM  is  assumed  to  be  completely 
oxidized  at  the  cathode  side.  The  methanol  mass  flux  is  constant  in 


MEM. 


the  ADL  and  MEM  due  to  no  electrochemical  reaction,  but  decreases 
in  the  ACL  as  a  result  of  methanol  consumption.  In  addition,  a  slight 
decrease  in  the  methanol  concentration  along  the  anode  flow  chan¬ 
nel  could  be  also  attributed  to  the  methanol  consumption. 

Oxygen  supplied  by  the  ambient  air  diffuses  through  the  CDL 
to  the  CCL  at  the  cathode.  The  oxygen  concentration  distribution 
in  the  cathode  porous  region  is  also  predicted  as  shown  in  Fig.  4. 
It  can  be  seen  that  the  oxygen  concentration  decreases  from  the 
CDL  to  the  CCL,  and  then  goes  down  non-linearly  across  the  CCL. 
Similarly,  the  oxygen  mass  flux  is  constant  in  the  CDL,  while  it 
decreases  in  the  CCL  due  to  the  oxygen  reduction  reaction.  In  con¬ 
trast  to  the  anode,  a  slight  increase  in  the  oxygen  concentration 
along  the  y-direction  is  observed.  And  the  main  reason  for  this 
slight  increase  may  be  considered  as  oxidation  of  the  permeated 
methanol  by  oxygen  in  the  CCL.  The  decrease  in  the  methanol  con¬ 
centration  along  the  y-direction  at  the  anode  results  in  a  decrease 
in  the  methanol  crossover  flux  leading  to  lower  oxygen  consump¬ 
tion  rate  at  the  cathode.  Thus,  the  oxygen  concentration  at  the  CDL 
surface  increases  according  to  the  boundary  condition  at  x  =  xcdl. 

3.3.  Pressure  distribution 

The  distribution  of  liquid  pressure  in  the  anode  porous  region  is 
shown  in  Fig.  5  where  the  liquid  pressure  near  the  interface  x  =  xadi 
is  magnified.  The  methanol  solution  is  fed  into  the  AFC  and  trans¬ 
fers  through  the  ADL  into  the  ACL  by  convection  and  diffusion.  It 
can  be  seen  from  Fig.  5  that  the  liquid  pressure  at  x  =  0,  that  is  the 
pressure  of  the  methanol  solution  in  the  AFC,  decreases  obviously 
along  the  flow  direction.  This  pressure  drop  results  from  the  friction 
resistance,  gravitation,  and  accelerated  flow  due  to  the  appearance 
of  C02  gas.  For  the  liquid  flow  in  y-direction,  the  flow  in  the  AFC 
dominates  the  flow  in  the  porous  region  due  to  the  assumption  of 
neglecting  the  flow  and  mass  transport  in  this  region.  Thus,  the  liq- 
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pressure  drop  is  about  700  Pa  from  the  ACL  to  the  ADL.  It  can  be 
expected  that  decreasing  liquid  pressure  in  y-direction  dominates 
the  decrease  in  the  gas  pressure  along  the  flow  direction,  while  the 
increased  capillary  pressure  resulted  from  the  decreased  liquid  sat¬ 
uration  (cf.  Fig.  8b)  determines  the  non-linear  increase  in  the  gas 
pressure  in  x-direction.  It  should  be  noted  that  the  liquid  and  the 
gas  pressure  gradients  are  different  at  the  interface  between  the 
ADL  and  ACL,  though  the  flow  is  continuous  across  the  interface. 
This  can  be  understood  from  Eqs.  (11)  and  (12)  considering  differ¬ 
ent  permeabilities  {I<  and  krg)  of  these  two  layers.  Furthermore,  it 
should  be  pointed  out  that  the  gas  pressure  is  lower  than  the  liq¬ 
uid  pressure  in  the  anode  porous  region  because  of  the  negative 
capillary  pressure  resulting  from  its  hydrophobicity. 

Water  is  produced  in  the  CCL  by  oxygen  reduction  reaction  and 
expelled  through  the  CDL  by  the  capillary  pressure  in  the  cathode. 
Fig.  7  shows  the  liquid  pressure  distribution  in  the  cathode  porous 
region.  The  cathode  liquid  pressure  is  constant  along  the  y-direction 
resulting  from  a  constant  cathode  gas  pressure  and  a  constant  liquid 
saturation  at  the  CDL  surface  facing  the  atmosphere.  However,  the 
calculation  shows  a  significant  drop  (1449  Pa)  in  the  liquid  pressure 
across  the  CCL,  which  is  mainly  resulted  from  a  very  low  permeabil¬ 
ity.  Then  the  liquid  pressure  decreases  linearly  through  the  CDL. 
Similar  to  the  anode,  different  liquid  pressure  gradients  appear  at 
the  interface  between  the  CDL  and  CCL  due  to  different  permeabil¬ 
ities,  and  the  liquid  pressure  is  also  higher  than  the  gas  pressure  (a 
constant  value  of  101300  Pa)  resulting  from  its  hydrophobicity. 


uid  pressure  in  the  ADL  and  ACL  decreases  along  the  flow  direction 
similar  to  that  in  the  AFC.  However,  the  liquid  pressure  along  the  x- 
direction  is  almost  constant  in  the  ADL,  and  decreases  slightly  in  the 
ACL.  The  slight  variation  of  the  liquid  pressure  in  the  porous  region 
may  be  attributed  to  the  positive  pressure  gradient  from  the  anode 
to  the  cathode  established  in  the  operation  cell  which  hinders  the 
liquid  convection.  Hence  diffusion  is  the  dominator  in  the  anode 
porous  region  resulting  in  the  slight  change  of  the  liquid  pressure. 
Following  the  methanol  arriving  the  catalyst  sites,  C02  gas  is  pro¬ 
duced  by  the  methanol  oxidation  reaction  in  the  ACL  and  transfers 
through  the  ADL  into  the  AFC  by  diffusion  and  convection.  The  gas 
pressure,  which  is  determined  by  the  liquid  pressure  and  capillary 
pressure,  cf.  Eq.  (13),  is  predicted  and  shown  in  Fig.  6.  The  gas  pres¬ 
sure  decreases  along  the  y-direction  and  varies  significantly  and 
non-linearly  across  the  porous  region.  The  calculation  shows  the 


3.4.  Liquid  saturation  and  relative  humidity  distribution 

The  predicted  liquid  saturations  in  the  AFC  and  in  the  anode 
porous  region  at  different  locations  in  the  y-direction  are  shown 
in  Fig.  8.  It  is  clear  that  the  liquid  saturation  in  the  AFC  decreases 
linearly  along  the  flow  direction  due  to  the  methanol  solution 
consumption  and  the  produced  C02  gas  accumulation.  The  liquid 
saturation  in  the  anode  porous  layers  keeps  almost  same  along  the 
y-direction,  which  is  mainly  caused  by  the  assumption  of  neglect¬ 
ing  flow  and  mass  transports  along  the  y-direction.  Furthermore, 
the  liquid  saturation  decreases  slowly  from  x  =  0  in  the  ADL  as  a 
result  of  one-way  transport  of  the  methanol  solution.  Then,  the  liq¬ 
uid  saturation  drops  dramatically  from  0.92  to  0.41  at  the  interface 
between  the  ADL  and  ACL,  and  decreases  slightly  in  the  ACL  due  to  a 
rather  small  capillary  pressure  gradient  (cf.  Figs.  5  and  6).  The  main 


Fig.  5.  The  distribution  of  liquid  pressure  in  the  anode  porous  region. 
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Fig.  6.  The  distribution  of  gas  pressure  in  the  anode  porous  region. 


cause  which  results  in  the  significant  drop  in  the  liquid  saturation 
at  the  interface  is  the  different  porosity  and  permeability  at  these 
two  layers  under  the  continuous  capillary  pressure,  especially  the 
permeability  of  the  ADL  being  about  two  orders  of  magnitude  larger 
than  that  of  the  ACL. 

Distributions  of  the  liquid  saturation  in  the  cathode  porous 
region  at  different  locations  in  the  y-direction  are  shown  in  Fig.  9a. 
The  cathode  liquid  saturation  almost  remains  same  along  the 
y-direction  due  to  the  constant  liquid  pressure,  cf.  Fig.  7.  Accom¬ 
panying  with  the  process  of  liquid  water  transferring  from  the  CCL 
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Fig.  7.  The  distribution  of  liquid  pressure  in  the  cathode  porous  region. 
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Fig.  8.  Distributions  of  (a)  liquid  saturation  in  the  AFC  and  (b)  liquid  saturation  in 
the  anode  porous  region  at  different  locations  in  the  y-direction. 


D.  Ye  et  al.  /  Journal  of  Power  Sources  192  (2009)  502-514 


511 


Fig.  9.  Distributions  of  liquid  saturation  in  the  cathode  porous  region:  (a)  at  different 
locations  in  they-direction;  (b)  at  y  =  0.0  mm  for  different  current  densities. 

to  the  outer  surface  of  CDL,  the  cathode  liquid  saturation  decreases 
non-linearly  across  the  CCL,  then  rises  from  0.035  to  0.058  at  the 
interface  between  the  CCL  and  CDL,  and  finally  decreases  linearly 
across  the  CDL.  The  jump  of  the  liquid  saturation  at  the  interface 
x  =  xad  can  also  be  attributed  to  the  continuous  capillary  pressure 
but  different  permeability  and  porosity  at  these  two  layers.  Fur¬ 
thermore,  the  liquid  saturation  distributions  through  the  cathode 
at  the  same  location  for  different  current  densities  are  shown  in 
Fig.  9b.  With  the  increase  in  the  current  density,  the  liquid  satura¬ 
tion  increases  in  both  the  CCL  and  CDL  due  to  the  increased  water 
resulting  from  the  oxygen  reduction  reaction  and  water  crossover 
from  the  anode.  Flowever,  it  is  noted  that  the  increase  in  the  liquid 
saturation  at  the  CCL  is  higher  than  that  at  the  CDL. 

The  evaporation  of  the  liquid  water  is  taken  into  account  in 
the  present  model.  The  relative  humidity  in  the  cathode  porous 
region  can  be  obtained  from  the  concentration  of  water  vapor  cou¬ 
pled  with  the  calculations  of  current  density,  oxygen  concentration, 
liquid  saturation  and  pressure.  Fig.  10  shows  distributions  of  the 
relative  humidity  in  the  cathode  at  current  densities  of  93.30  A  m-2, 
364.61  Am-2,  and  720.78  Am-2,  respectively.  The  relative  humid¬ 
ity  decreases  across  the  cathode  porous  region  from  the  CCL  to  the 
CDL  for  all  current  densities.  The  relative  humidity  increases  with 
increasing  current  density  comparing  Fig.  lOa-c.  It  can  be  expected 
that  a  large  current  density  leads  to  a  high  production  rate  of  the 
liquid  water  and  a  large  water  crossover  flux,  hence  a  high  interfa¬ 
cial  transfer  rate  between  liquid  water  and  vapor.  In  addition,  the 
relative  humidity  along  y-direction  is  almost  uniform  as  the  result 


of  the  constant  liquid  saturation  in  y-direction,  which  can  be  easily 
understood  with  the  help  of  Fig.  9a. 

3.5.  Methanol  crossover 

It  is  known  that  methanol  crossover  will  definitely  deteriorate 
the  cell  performance.  As  mentioned  before,  diffusion,  electro¬ 
osmosis  and  hydraulic  pressure  difference  determine  the  methanol 
crossover  flux.  Fig.  11  shows  the  variation  in  the  methanol  crossover 
flux  with  the  current  density  for  various  methanol  concentrations. 
For  a  given  current  density,  the  methanol  crossover  flux  increases 
with  increasing  methanol  concentration.  For  a  given  methanol  con¬ 
centration,  it  increases  with  increasing  current  density  for  the 
concentrations  of  2  M  and  4  M,  while  it  increases  slightly  at  low  cur¬ 
rent  densities  and  decreases  at  high  current  densities  for  1  M.  For 
further  understanding  the  mechanism  of  the  methanol  crossover, 
the  total  methanol  crossover  flux  and  its  three  contributors  at  dif¬ 
ferent  feed  methanol  concentrations  are  predicted  by  the  present 
model  as  shown  in  Fig.  12.  It  is  can  be  seen  that  for  both  4M 
and  1  M  feed  concentrations,  diffusion  predominates  the  methanol 
crossover  at  low  current  densities,  while  electro-osmosis  is  the 
dominator  at  high  current  densities.  The  weakening  effect  of  the  dif¬ 
fusion  on  the  methanol  crossover  flux  can  be  explained  by  the  fact 
that  a  high  current  density  leads  to  a  high  methanol  consumption 
rate  hence  a  low  methanol  concentration  at  the  interface  between 
the  ACL  and  MEM  for  either  methanol  concentration.  Flowever,  dif¬ 
ferent  enhancing  effect  of  the  electro-osmosis  on  the  methanol 
crossover  at  various  methanol  concentrations  can  be  seen  from 
Fig.  12.  The  methanol  crossover  flux  driven  by  the  electro-osmosis 
significantly  increases  with  increasing  current  density  at  4  M,  while 
it  first  increases  and  then  decreases  at  1  M.  The  variation  in  the 
methanol  crossover  flux  driven  by  the  electro-osmosis  is  due  to  the 
combined  effect  of  the  increased  current  density  and  the  decreased 
methanol  concentration  at  the  interface  x  =  xad.  In  addition,  the 
hydraulic  pressure  difference  has  a  negative  effect  on  the  methanol 
crossover  due  to  higher  liquid  pressure  at  the  cathode,  but  the  vari¬ 
ation  of  the  absolute  value  in  the  methanol  crossover  flux  driven  by 
the  pressure  difference  is  similar  to  that  by  electro-osmosis  for  both 
4  M  and  1  M.  This  result  is  also  attributed  to  a  combined  effect  of 
the  increased  cathode  liquid  pressure  and  the  decreased  methanol 
concentration  at  x  =  xacl.  It  can  be  concluded  that  decreased  con¬ 
centration  at  the  interface  x  =  xad  plays  a  critical  role  in  methanol 
crossover  for  low  feed  methanol  concentrations,  which  leads  to  the 
decrease  in  the  methanol  crossover  flux  at  high  current  densities. 

Taking  account  of  the  methanol  crossover,  the  fuel  utilization 
efficiency  in  the  DMFC  can  be  calculated  by: 

=  /TV  (34) 

Fig.  13  shows  the  fuel  utilization  efficiency  for  various  methanol 
concentrations.  The  fuel  utilization  efficiency  increases  with 
increasing  current  density  because  of  a  decreased  ratio  of  the 
parasitic  current  density  to  the  current  density,  even  though  the 
methanol  crossover  flux  increases  simultaneously  at  high  feed 
methanol  concentrations.  In  addition,  high  fuel  utilization  effi¬ 
ciency  at  low  methanol  concentrations  is  observed  due  to  small 
methanol  crossover  flux  as  discussed  before. 


3.6.  Water  crossover 


The  net  water  transport  coefficient  is  widely  used  to  describe 
the  water  crossover  flux  through  the  MEM  [42,43].  According  to  Eq. 
(18),  the  net  water  transport  coefficient,  aw,  can  be  obtained  by: 


aw  =N™em  j  =  nd  w 


F  A  /f(Pl,c-Pl,a) 
/  Mw  /-'Tmem 


(35) 
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Fig.  10.  Distributions  of  relative  humidity  in  the  cathode  at  different  current  densities:  (a)  93.30  A  rrr2;  (b)  364.61  Am-2;  (c)  720.78  A  irr2. 


The  electro-osmotic  drag  coefficient  of  water,  ndw,  is  a  constant 
at  a  given  temperature  (see  Table  4).  Thus,  the  net  water  trans¬ 
port  coefficient  depends  on  the  current  density  and  the  hydraulic 
pressure  difference  across  the  MEM.  A  positive  aw  represents  a  net 
water  transport  from  the  anode  to  the  cathode,  while  a  negative 
aw  denotes  a  reversed  transport.  Fig.  14  shows  the  predicted  net 
water  transport  coefficient  and  its  contributors  for  4  M  feed  con¬ 
centration.  It  is  can  be  seen  that  the  net  water  transport  coefficient 


increases  rapidly  at  low  current  densities,  and  then  tends  to  a  mild 
increase  for  current  densities  higher  than  200  Am-2.  This  can  be 
understood  by  the  variation  in  the  water  crossover  flux  driven  by 
electro-osmosis  and  hydraulic  pressure  difference.  It  is  clear  that 
the  water  crossover  flux  is  dominated  by  the  electro-osmosis  and 
the  hydraulic  pressure  difference  across  the  MEM  at  low  current 
densities,  while  electro-osmosis  plays  an  essential  role  at  high  cur¬ 
rent  densities. 
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Fig.  11.  The  variation  in  the  methanol  crossover  flux  with  the  current  density  for 
various  methanol  concentrations. 


Fig.  13.  The  fuel  utilization  efficiency  for  various  methanol  concentrations. 


In  addition  to  the  water  crossover  flux,  the  total  water  flux  at  the 
cathode  is  predicted  as  follows: 

N“tal  =  ACem  +  N°RR  +  N°rpm  ,  (36) 

where  N™em  denotes  the  water  crossover  flux  from  the  anode  to 
the  cathode,  N°RPM  represents  the  water  flux  produced  by  the  ORR 


Current  density  /  A  m'2 


in  the  cathode  given  by: 

N°RR=^=,  (37) 

and  N°RPM  is  the  water  flux  generated  by  the  oxidation  reaction  of 
permeated  methanol  (ORPM)  calculated  from: 

N°RPM  =  lf-r  (38) 

The  variation  in  the  total  water  flux  and  its  contributors  with  the 
current  density  at  4M  feed  concentration  is  shown  in  Fig.  15.  It  is 
noted  that  a  significant  increase  in  the  total  water  flux  is  observed 
with  increasing  current  density,  which  can  be  understood  from  the 
characters  of  its  three  contributors.  First,  the  water  crossover  flux 
increases  evidently  with  the  increase  in  the  current  density  as  men¬ 
tioned  above.  Second,  the  water  flux  produced  by  the  ORR  is  directly 
proportional  to  the  current  density  according  to  Eq.  (37).  At  last,  the 
water  flux  from  the  ORPM  is  in  direct  proportion  to  the  parasitic  cur¬ 
rent  density  which  increases  with  increasing  current  density  (cf. 
Fig.  11).  Moreover,  it  can  be  concluded  from  Fig.  15  that  the  total 
water  flux  at  the  cathode  is  originated  primarily  from  the  water 
generated  by  the  oxidation  reaction  of  the  permeated  methanol  at 
low  current  densities.  Flowever,  the  water  crossover  flux  is  the  main 
source  of  the  total  water  flux  at  the  cathode  once  the  current  density 
is  higher  than  300  Am-2. 


Fig.  12.  The  total  methanol  crossover  flux  and  its  three  contributors  at  different  feed  Fig.  14.  The  net  water  transport  coefficient  and  its  contributors  at  4  M  feed  concen- 
methanol  concentrations:  (a)  4  M;  (b)  1  M.  tration. 
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Fig.  15.  The  variation  in  the  total  water  flux  and  its  contributors  with  current  density 
at  4  M  feed  concentration. 

4.  Conclusions 

A  two-dimensional  two-phase  mass  transport  model  is  devel¬ 
oped  for  a  vertical  air-breathing  direct  methanol  fuel  cell,  where 
the  discontinuity  of  liquid  saturation  at  the  interface  between  the 
diffusion  layer  and  catalyst  layer  as  well  as  the  mass  transport  in  the 
catalyst  layer  is  considered.  The  predicted  cell  performance  agrees 
well  with  the  experimental  data  within  a  wide  current  density 
range.  The  typical  two-phase  flow  and  mass  transport  distributions 
at  both  the  anode  and  the  cathode,  the  methanol  crossover,  the  net 
water  transport  coefficient,  the  water  crossover  flux,  and  the  total 
water  flux  at  the  cathode  and  their  contributors  are  predicted.  The 
main  results  of  the  present  model  are  summarized  below: 

(1 )  The  methanol  concentration  decreases  across  the  ADL,  ACL  and 
MEM  and  goes  down  slightly  along  the  flow  direction  in  the 
AFC,  and  the  oxygen  concentration  decreases  from  the  CDL  to 
CCL  but  increases  along  the  flow  direction  of  methanol  solution. 

(2)  The  anode  gas  pressure  increases  greatly  across  the  ADL,  while 
the  cathode  liquid  pressure  drops  significantly  across  the  CCL. 

(3)  The  liquid  saturation  decreases  significantly  across  the  ADL  and 
the  CCL,  and  it  drops  to  a  smaller  value  from  the  diffusion  layer 
to  the  catalyst  layer  at  both  the  anode  and  cathode. 

(4)  The  weakening  effect  of  diffusion  and  the  different  enhancing 
effect  of  electro-osmosis  on  the  methanol  crossover  are  found 
for  various  methanol  concentrations.  Diffusion  dominates  the 
methanol  crossover  at  low  current  densities,  while  the  electro¬ 
osmosis  is  the  dominator  at  high  current  densities. 

(5)  Water  transport  through  the  MEM  depends  on  electro-osmosis 
and  hydraulic  pressure  difference  across  the  MEM  at  low  cur¬ 
rent  densities,  while  electro-osmosis  plays  a  critical  role  in  the 
water  crossover  at  high  current  densities.  The  total  water  flux  at 
the  cathode  is  originated  primarily  from  the  water  generated  by 


the  oxidation  reaction  of  the  permeated  methanol  at  low  cur¬ 
rent  densities,  while  the  water  crossover  flux  is  the  main  source 
of  the  total  water  flux  at  high  current  densities. 
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